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Abstract 

We extend Gegenbauer Polynomials technique to evaluate a class of complicated Feynman 
diagrams. New results in the form of si^-hypergeometrical series of unit argument, are pre- 
sented. As a by-product, we present a new transformation rule for 3 F 2 -hypergeometric series 
with argument —1. 
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Though Quantum Field Theory already has a long history and a number of different ap- 
proaches, Feynman diagrams (FD) are still the main source of its dynamical information. The 
necessity to know more exactly about the characteristics of physical processes and specific 
QFT functions themselves stimulates the calculation of radiative corrections of ever higher or- 
der. This, in its turn, leads to new methods of calculation (see, for review, Ref.0). Fifteen 
years ago the method based on the expansion of propagators in Gegenbauer series (see has 
been introduced in []3], f|]. One has shown that by this method the analytical evaluation of 
counterterms in the minimal subtraction scheme at the 5-loop level in any model and for any 
composite operator was indeed possible. The Gegenbauer Polynomial (GP) technique has been 
applied successfully for propagator-type FD in many calculations (see jl]- ||). In the present 
paper we present a new development of the GP technique and apply it to the evaluation of a 
class of master two-loop diagrams containing the vertex with two propagators having the powers 
1 or D/2 — 1 in D- dimensional space. 

Thoughout the paper we use the following notation. The use of dimensional regularization is 
assumed. All the calculations are performed in the coordinate spaceQ of dimension D = 4 — 2e. 
Any diagram is converted from p-space to x-space by applying the Fourier transform (marked 

by the symbol =^>) 

d D pexp(ipx) d/2 d-2u a o(") / x T(D/2-u + n) 



(p 2 )" J (p 2 )" (x*)D/2-*> ^ > Y{y) 

(r is the Euler T- function) or by considering the dual diagram (see, for example, The 
trasformation to the dual diagram will be denoted by =. 

The dual diagram is obtained from the original one by replacing p, by Xj into corresponding 
integral and using the rules of correspondence between the graph and the integral, as in the 
x-space. Thus, after applying the Fourier transform we obtain a different D-space integral 
but the x-space FD similar graphically to corresponding p-space one. Considering the dual 
transformation, we do not change the integrand, except for the replacement Pi — > Xj, but 
obtain another graphical representation. Using of the both transformations leads to define the 
number of really independent FD (see section 2). 

Note that contrary to || we analyse FD directly in x-space which allows us to avoid the 
appearence of Bessel functions. It is possible in the case of propagator-type FD, because their 
dependence on a single external momentum (in p-space) or on a single external coordinate is 
power-like and known beforehand. The point of interest is the coefficient function Cf, which 
depends on D = 4 — 2e and is a Laurent series in e. If some diagram has Z-lines^ with the 
corresponding indices a{, i = {1, ...,/} and L-loops in p-space, then its coefficient functions in 
p-space and x-space are connected by 

Cf(p — space) = — — °^ - — — • Cf(x — space) and Ct(p — space) = Ctix — space) 

i - DI. 2) 

in the case of the Fourier and dual transform, respectively. 



2 Note that the choice of x-space is not crucial. If the reader prefers to work in p-space he can consider the 
variable {xi, yi, ...} as some momenta in p-space. 

3 With the lines of the graphs there are associated power-law factors of the type (x 2 )~ a , in which a is called 
the line index. 
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Following [(J we introduce the traceless product (TP) i' 11 '"' 1 " connected with the usual 
product x Ml ...x Mn by the following equations (see 0, |fl) 



q V "-■( l) p r(n p+X) „ lft2 u 2p _ lfl2p 2 P M2 P+ i ™Mn 
*£ 2*p!(n-2p)!r(n + A) 9 ~' 9 X * - X 



nWfn — In 4- A 4- 1 "I 
^ (2) 2 Pp!(n - 2p)!r(n - p + A + 1) 

where hereafter A = D/2 — 1. 

Comparing Eqs. (|A.4j ) and ([!]), we obtain the following relation between TP and GP 



^l.-Mn ^.-Mn = ^ !r ( A ) C*(XZ) (x'z'f 2 (2) 

2 n r(n + A) nV J y ' y ' 

We give also the simple but quite useful conditions: 

Z w-Mn x m—iMi _ 2 Mn x^ 1 '"^" = z fll '"^ n x >11 x^ n (3) 

which follow immediately from the TP definition: x^i-m-m-v™ = o. 

The use of the TP ^Mi-Mn makes it possible to ignore terms of the type g^ 1 that arise upon 
integration: they can be readily recovered from the general structure of the TP. Therefore, 
in the process of integration it is only necessary to follow the coefficient of the leading term 
x Ml ...x^ n (see also section 3). The rules to integrate FD containing TP can be found, for example 

lir) D ): 



in || 0. For a chain we have (hereafter Dx = (d D x)/(2 

T fll...fJ.n 1 „.^\...^ln 

X 2a {x-y) 2 P (4tt) D / 2 y2( a +/3-A-l) ^ 



where 



a n+m (a + /3- A - 1) 

Note that in our analysis it is necessary to consider more complicate cases of integration, 
when the integrand contains G functions. Indeed, using the Eqs.( [A.2|) and (fj), we can represent 
the propagator (x\ — X2)~ 2X into the following form: 



E 



2 n r(n + A) 



Lll-IMl 111... fin 
th 1 Ju • 



^^0(x 2 2 -xl) + (x\< — >xf) 



(5) 



Using the GP properties from Appendix A and the connection (0) between GP and TP, we 
obtain the rules for calculating FD with the B-terms and TP. 

1. The rules have the following form: 

x ^ n ^ 1 ^ B(m,n\(3,\) 

X trial „ _ ,A2/3^ X y > (A^\D/2 „2(a+/3-A-l) 2^ 



x 2a (x-y) 2 P v tf/ (4tt) d / 2 y^+P-^ 1 ) ^ Q m + a + f3-\-l 

(/3 =A) 1 y^ 1 -^ _L 1 

(Att) d / 2 7^ r(A) (a - l)(n + A) 



(6) 



2 



Dx 



(J3=X) 



X 



x 2a (x - y) 2 P 



e(y 2 -x 2 



1 



y 



fll. ..fin 



(47r) D/2 y 2(a+P-X-l) E q 
1 



B{m,n\(3,X) 
m + n — a + A + 1 



(Att) d / 2 y 2 ^- 1 ) T(X)(n + X + l-a)(n + X) 



(7) 



where 



B(m,n\fi, X) 



T(m + fi + ri) 



T(m + P-X) 



m\T(m + n + l + X)T(J3) T(J3 - A) 

Note that the x-space index of an ordinary line is A, where all our formulae are simplified. 
Hereafter we add this specific case (3 — X to our rules. 

The sum of above diagrams does not contain 9-terms and should reproduce Eq.(f|). To 
compare the r.h.s. of Eqs.(P,0) and the r.h.s. of Eq.([|) we use the transformation of 3_F 2 - 
hypergeometric function with unit argument (see |flO[|): 



3 F 2 (a,b, c;e,f; 1) 



r(i-o)r(e)r(/)r(c-6) 

T(e - 6)r(/ - 6)r(l + b - a)r(c) 
3^(6,6- e + 1,6- / + 1; 1 + 6- c, 1 + 6- a; 1) + (b 



-> c 



(8) 



When e = 6 + 1, the 3 .F 2 -function can be represented as the sum of another 3 F 2 -function 
and a term containing only T-function products: 



E 

k=0 



T(k + a)T(k + c) 1 
k\T{k + f) k + b 



T(l - a)r(a) 



T(f-c)T(l + c-f) 



E 

fc=0 



r(o)r(i -a)r(6)r(c-6) 

r(/-6)r(l + 6-a) 

r(k + c- f + i)r(fc + c ) i 

k\T(k + 1 + c - a) k + c-b 



(9) 



This is the case (when k = m, b = a + (3 — A — 1, c = n + /?) to compare Eq.(f|) and the sum of 
eqs.(§0). 

Analogously to Eqs.(^) and (0) we have more complicate cases 4 : 



/ 



X' 



HX—Hn 



x 2a (x - y) 2 ? 



6(x 2 



(Att) d / 2 



y 



H\. ../Mi 



Q^ 2 -^) A n,0 (a g) 
2(o+/9-A-l) A 



m=0 
G9=A) 



^(a+jS-A-l) 
1 1 



2 m+P+n 



B(y 2 - + 



+ 



(4tt) d / 2 T(A) y 
1 1 /6(z 2 - y : 



m + a + /5 — A — 1 v y z 
1 9(y 2 -z 2 ) 

•2.n+A e(?/ 



m — a + n + l + A 



^(a-l) „ _|_ \ \ a _ 1 



- 2 - ^ 



n + 1 + A — a 



(10) 



X 



fll-.-fln 



x 2a (x - y) 2f3 



®{z 2 -y 2 ) 



(4tt)^/ 2 



111... fir, 



Q{z 2 -y 2 ) Anfi 



y 



2(a+/3-A-l) 



A"-°(a,/9) 



J The full set of rules will be presented in 
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~ B(m,n\P,\) (,y*_ 

2^ v 2(a+d-X-i\ \ \ z 2 



m=0 
(/3=A) 1 



1 



&(z 2 -y 2 ) 



y 



(4tt) d / 2 r(A) 

i i /e(z 2 



m + a + (3 — X — 1 
1 &(z 2 - 



2 m+f3+n 



2,2(0-1) ( a _ i)( n + x + 1 - a) 



z 2(a-l) n _|_ \ 



2 n+X 



a — 1 



n + 1 + A — a 



m — a + n + l + X, 



One can easily see that the sum of the above diagrams lead to results identical to ([|). 

2. The aim of this article is to study a class of master two- loop diagrams containing the 
vertex with two propagators having index 1 or A. 

Consider in the x-space the following general diagram 



/ 



DxDy 



y 2a (z — y) 2t (z — x) 2/3 x 2 " l (x — y) 



2 s 



= J(a,t,P,"/,s) 



and restrict ourselves to the FD A(a, j3, 7) = J (a, A, P, 7, A), which is the one of FD of interest 
for us here. It is easily shownP] that 



C f [A{a,P,i)\ = C f [A{a,<T,i)\ = C f [J(j, A, X,a,a)) = C>[J(<7,7, A, A, /?)], 



Doing Fourier transformation of both: the diagram A(a, P, 7) and its solution in the form 
Cf[A(a, P, 7)](z 2 ) _CT , where hereafter t — X + l—t, t = {a, j3, 7, a, ...}, and considering the new 
diagram as one in the x-space (i.e. making the dual transformation) we obtain the relation 



C f [A(a,P, 7 )) 



a Q (8) 



C f [J(a,l,P,%l)} 



ao(A)a (a)a (/3)ao( 7 )«o(o r ) Cf[J((3, 1,6,7, 1)] 



(13) 



between the diagram, which contains the vertex with two propagators having the index A, with 
the similar diagram containg the propagators with the index 1. 

Repeating the manipulations of subsection 1 we can obtain the following relations: 

C f [J0, 1, a, 7, 1)] = C f [J(j3, 1, a, 7, 1)] = C f [J(l, 1, 7, a, a)} = C f [J(a, 1, 1, 7, P)} (14) 

Thus, we have obtained the relations between all diagrams from the class introduced in the 
beginning of this section. Hence, it is necessary to find the solution for one of them. We prefer 
to analyse the diagram A(a,P,j), that is the content of the next section. 



3. We calculate the diagram A(a,P,j) by the following wayQ 

M! ^ A-L\n -+- A) f Z Hl-Vn yUl-Hn |-0(^ 2 

SiyOi, p, 7J 



£ , DxDy 



n=0 



nW(X) 



x 2l {z — x) 2f3 y 2a [x — y) 



2X 



y 2 ) Q{y 2 -z 2 ) 



z 



2(n+A) 



+ 



y 



2(n+X) 



5 For example, we can consider new FD which is the result of gluing of A(a, (3, 7) and the propagator with index 
ct + A, where a = 3 + A— (a + f3 + ^) and A — > 0. This new FD has a pole 1/A with the corresponding coefficient 
independent from the external moment. Hence, this coefficient does not depend from which propagator of new 
FD is external and, thus, the values of these different (nongluing) FD are connected. 

(n) 

6 The symbol = markes the fact that the equation (n) is used on this step. 
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<@P 1 J_ 1 ~ 2T(rc + A) 

(An) D / 2 r(A) a-l±i 



ra=0 



n!r(A) 



e(^ 2 - x 2 



+ 



e(x 2 -z 2 



Dx 
1 



2' 



A t i---A t n^A t i---Ai 



1 



&(z 2 — X 



\ + n + 1 — a 



+ 



e(x 2 



^ z 2{n+\) x 2{a-l) x 2(n+\) z 2(a-l) 1 \ _|_ n _|_ a _ \ ^2(n+A+a-l) x 2(n+A+a-l) ■ 

Evaluating the second integral and using the equation (see ( |A.1| ) and @) 

T(n + 2A)r(A) 



we have from Eqs.(|6|) and (|^): 

Cf[A(a,p n )) = 
1 



1 



2 n r(2A)r(ra + A) ' 



1 ^ T(n + 2A) ~ 



(4tt)^ r(A) a - 1 „f «!r(2A) m=0 

1 1 

+ 



A + n + 1 — a m + n + 2 — « — 7 + A m + n + 7 + /3 — 1' 
1,1 1 



+ 



X + n + a — 1 m + n + A + 1 — 7 m + n + a + 7 + (3 — 2 



(15) 



(16) 



The r.h.s. of Eq. (jig ) is already in an easy-to-use form, because the sums converge quite fastly 
and thus they can be easily evaluated numerically. However, the aim of our investigations is to 
obtain the simplest analytical expresion for the initial diagram. As we will see in this section 
it is 3 F 2 -hypergeometric functions with unit argument. 

To our purposes it is useful to transform the parts 

1 1 . . 1 1 



m + n + 2 



a 



7 + A J> ~m + n + 'y + (3 — 1 



'm + a + j + [3 — A 



+ 



m + 1 — 7' 



1 



+ 



1 



1 



1 



m + n + A + 1 — 7 m + n + a + 7 + /5 — 2' 



l m + 7 + /3 — 1 — A m + 2 — a — 7' 



It is can be done by using of Eq.([5]) or by represent at ing Q(x 2 — x 2 ) terms into FD of the 
r.h.s. of Eq.(|15D as 1 — Q(x 2 — xf). Then the result for Cf[A(a, f3, 7)] transforms to 



C f [A(a,[3 n )} 



(4tt) d T(A) a - 1 



I — I 



where 



and 



E 



T(n + 2X) 



A + n + 



l — ( A n >°(a - 1 + 7, (3) + A n '°(n + A + 7, (3) 

■f 1 - « V 



n=o nir(2A) 

• (^'°(7, P) + A n >°(n + a + A + 7 - 1, /?)) 



A + n + 



~ £(n + 2A) ~ 



- nT(2A) 



m=0 



1 



1 



1 



A + n + 1 — a \m + a + (3 + 7 — A — 2 
1 1 



m + 1 — 7/ A + n + a — 1 \m + /3 + 7 — A — 1 m + 2 — a — 7 



(17) 
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Apply Eq. (^) with k = n, 6 = A + 1 — a and with k = m, b = a + (3 + r y — A — 2 and 
6 = 1 — 7 in the r.h.s. of Eq.(^). After some algebra the double-sum terms cancel one of 
another and we obtain the following form for the part / 



r(i - p)r{x + 1 - a)r(A - 1 + a)r(i -p + A)r(i - 7 )r(« + /3 + 7- A-2) 



- E 



T(2A)r(2 + A - a - /3)T(a + 7 - l)r(2 + A - 7 - p)T(a + (3 - A - 1) 
T(n + 2A) (-l)T(l-/?) 1 



, ): „ n!r(2A) r(/3-A) A + n + a-1 
r(a + /3 + 7 - A - 2)r(2 - a - 7) 



(19) 



+ 



r(i- 7 )r(/? + 7 -A-i) 



r(3 - a - p - 7 - n)T(a + 7 + A - 1 + n) T( 7 - A - n)T{2 + 2A - /3 - 7 + n) 



Starting from the term ~ (n + A — 1 



,-1 



in the r.h.s. of Eq.(|18D and repeating the above 



analysis, we can obtain another form of /: 

~ r(i - p)T{X + 1 - «)r(A - 1 + «)r(i -p + A)r(2 -a- j)T{p + 7 - A — 1) 



+ E 



T(2A)r(A + a - p)T(i)T(3 + A - a - 7 - (3)T{(3 - a - A + 1) 
T(n + 2A) (-l) n r(l-/5) 1 



n f n!r(2A) r(^-A) 

r(/3 + 7 -A-i)r(i- 7 ) 



A + n — a + 1 

r(2 - a - 7)r(a + /3 + 7- A-2) 



(20) 



r(2 - /3 - 7 - n)r(7 + A + n) T(a + 7 - 1 - A - n)r(3 + 2A-a-/5-7 + n) 



Applying Eq.(^) to one of the two terms contained in the sum in the r.h.s. of Eqs. (|19D or 
(PD|) we can have also two additional representations for /: 



T(l - p)T(X + 1 - P)T(X + 1 - a)r(A - 1 + a)r(l - 7)r(2 - a - 7) 



r(2A)r(a)r(i - a)r(2 + a - « - /3)r(3 + a- «- 7 -/?) 



+ E 



r(n + 2A) (-l)«r(l-/3) 



n^ ^(2A) 

1 



r(a + /3 + 7 - A - 2)r(2 - a - 7) 
r(3 + 2A - a - (3 - 7 + n)Y{a + 7 - A - 1 - n) 



A + n — a + 1 



rog-A) 

r(i- 7 )r(/? + 7- a - 1) 1 

r(7 - A - n)r(2 + 2A-/5-7 + n)'A + n + a 



(21) 



r(i - /3)r(A + 1 - /3)r(A + 1 - «)r(A - 1 + a)r(/3 + 7 - a - i)r(a + 7 + /3-A-2) 



r(2A)r(a)r(i - a)r( 7 )r(A + 7-1) 



r(n + 2A) {-l) n T{l-l3) 



+ n5 n\T{2\) 



r(p - A) 



1 



r(/? + 7-A-i)r(i- 7 ) 

r(A + 7 + n)r(2 - /3-7-n)"A + n- a + l 



r(2 - a - 7)r(a + /3 + 7- A-2) 
r(3 - a - 7 - P - n)T(a + A + 7 - 1 - 



n 



A + n + a — 1 



(22) 



Eqs. (|l9|)-(p^) are equal to one another because they are connected by the transformation 
(§). For concrete values of a, P and 7 the more convenient one of them may be used. 

Thus, a quite simple solution for A(a,P, 7) is obtained^]. In next section we will consider 
the important special case of these results. 



7 Before our studies, the possibility to represent Cf[A(a, (3, 7)] as a combination of 3-F 2 -hypergeometric 
functions with unit argument, has been observed in B. 
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4. As a simple but important example to apply these results we consider the diagram 
J(l, 1, 1, 1, a). It arises in the framework of a number of calculations (see JT1J - JT5|] ) . Its 
coefficient function 1(a) = C/[J(1, 1, 1, 1, a)} can be found as follows 

1(a) I ^(iKW c f [J(\, A, A, A, a)] and C f [J(X, A, A, A, a)] = C f [A(a, 3 - A - d, A)] 
ao(a + 2 — 2A) 

The latter equation may be obtained by analogy with (|Hp. 
From Eqs. (|17D and (^) we obtain 



2 r 2 (A)r(A-a)r(q + l-2A) 
W (4tt)^ r(2A)r(3A - a - 1) 1 ' 

r 2 (l/2)r(3A-a-l)r(2A-a)r(a + l-2A) + g Y(n + 2A) 



T(A)r(2A + 1/2 - a)r(l/2 - 2A + a) ^ T(n + a + 1) n + 1 - A + a 

Note that in [|11|] Kazakov has got another result for 1(a): 

2 r 2 (A)r(l-A)r(A-a)r(q + l-2A) 
1 J " (4vr)^ r(2A)r(a)r(3A - a - 1) 1 J 

T(A)r(2 - A)r(a)r(3A - a - 1) ~ T(n + 2A) / 1 



T(2A-l)r(3-2A) y r(n + 2-A) \n + l-X + a n + 2A - a 

From Eqs. (|23| ) and ( p4[ ) we obtain the transformation rule for 3F2-hypergeometric function 
with argument —1: 

, , , s , T(2-a)T(b + a-l)T(b-a)T(l + a-b) , . 

- - 1 ~ a - ■ 3^2(20, 6, 1; 6 + 1, b + a; 1) - — — ■ 3^2(20, 1 + a - 6, 1; 2 + a - 6, 2 - a; -1), 

b + a — 1 1 + a — b 

where a = A and b = 1 — A + a are used. 



Equation (25) has been explicitly checked at a = 1 and b = 2 — a (i.e. A = 1 and a = 1), 
where the 3 F 2 -hypergeometric functions may be calculated exactly. We cannot directly prove 
Eq.(p5|) at arbitrary a and b values: the general proof seems to be non-trivial. Note that it is 



different from the equations of |T6| and may be considered as a new transformation rule. 



5. As a conclusion we note that we obtained the results for the class of master two-loop 
diagrams containing the vertex with two propagators, having indices 1 or \ = D/2 — 1. As a 
by-product of the above-mentioned studies, the new transformation rule of 3F2-hypergeometric 
series with argument —1, was obtained. It may be used together with other rules (see Ref . [|T6|j , 
for example). 

We would like to note also, that there is only one other diagram 

7) = J( a , A,/?, A, 7) 

with two indices A, which does not fall into the above-considered FD class. The diagram 
J(a, A, (3, A, 7) can be obtained from B(a,/3, 7) by a analysis similar to ([U|) and it cannot be 
considered as the independent one. The incorporation of the diagram B(a, (3, 7) in our analysis 
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invites futher investigations and it is the subject of the next publication 

The author is grateful to Prof. A.N.Vasil'ev for a stimulation of this research and Prof. 
D.Broadhurst for information about a possibility to represent FD A(a, /3, 7) as a combination 
of 3 F 2 -hypergeometric functions with unit argument (see |J). 
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Dr. A.L.Kataev, Prof. D.I. Kazakov, Dr. N.A.Kivel, Dr. A.S. Stepanenko and Prof. A.N.Vasil'ev 
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A Appendix 

Here we present useful formulae to use of Gegenbauer polynomials. Following |3|, |J, D-space 
integration can be represented in the form 

d D x = -SD-i(x 2 ) X dx 2 dx, 

where x = x/y/x*, and Sd-i = 27r A+1 /T(A + 1) is the surface of the unit hypersphere in Rd- 
The Gegenbauer polynomials C^{t) are defined as 0, ^[ 

(l-2rt + rV = E C n (ty (r<l), Cj(l) = ^t^f (A.1) 
whence the expansion for the propagator is: 

J2\ n / 2 



oc 



(A.2) 



n=0 

where 

'l, if y > 
0, if y < 

Orthogonality of the Gegenbauer polynomials is expressed by the equation (see 

A 



e(y) 



C^fifa) *a C^(f 2 f 3 ) = ~—5n C^x\x 3 ), (A.S) 

n + A 



where 8™ is the Kroneker symbol. 

The following formulae are useful (see || ||): 

CniX) ~ h (n-2p)\p\T(6) 

^ = E^ g ) (A . 4) 

n! p>o p!r(n - p + d + 1) 

Substituding the latter equation from ( |A.4j ) for 5 = A to the first one, we have the following 
equation after the separate analysis at odd and even n: 



[n/2] 



(n - 2k + A)T(A) T(n + 5 - k)T(k + 5 - A) 



n ~ fe = p fc!r W r(n + A + i-fc)r(5-A) 1 bj 

This equation is an example of the equations used in our analysis. Note that alternatively 
to ( |A.5| ) we can use the expansion of (x\ — a^) -25 similar to (|A.3|) but directly in a series of GP 
C^Xi-xi) (seeffl). 
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